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ABSTRACT 

The “4-day wave” is an eastward moving quasi-nondispersive feature with period near 4 days occurring near 
the winter polar stratopause. This paper presents evidence of the 4-day feature in Microwave Limb Sounder 
(MLS) temperature, geopotential height, and ozone data from the late southern winters of 1992 and 1993. Space- 
time spectral analyses reveal a double-peaked temperature structure consisting of one peak near the stratopause 
and another in the lower mesosphere, with an out-of-phase relationship between the two peaks. This double- 
peaked structure is reminiscent of recent three-dimensional barotropic/baroclinic instability model predictions 
and is observed here for the first time. The height variation of the 4-day ozone signal is shown to compare well 
with a linear advective-photochemical tracer model. Negative regions of quasigeostrophic potential vorticily 
(PV) gradient and positive Eiiassen-Palm flux divergence are shown to occur, consistent with instability dynamics 
playing a role in wave forcing. Spectral analyses of PV derived from MLS geopotential height fields reveal a 
4-day signal peaking near the polar stratopause. The three-dimensional structure of the 4-day wave resembles 
the potential vorticily “charge” concept, wherein a PV anomaly in the atmosphere (analogous to an electrical 
charge in a dielectric material) induces a geopotential field, a vertically oriented temperature dipole, and cir- 
culation about the vertical axis. 



1. Introduction 

The 4-day wave is a ubiquitous feature in the polar 
winter upper stratosphere. It was first observed in tem- 
perature data as a strong 4-day signal in zonal wave- 
number 1 (Venne and Stanford 1979). Later studies 
showed that it consists of waves 1 through at least 4 all 
moving with the same phase speed, such that the period 
of wave 1 is near 4 days. Synoptic plots of temperature 
reveal a quasi-nondispersive “warm pool” of air, which 
rotates eastward around the winter pole near the stra- 
topause with a 4-day rotation period (Prata 1984; Lait 
and Stanford 1988b; Lawrence et al. 1995). Numerous 
observational (Venne and Stanford 1979, 1982; Prata 
1984; Lait and Stanford 1988b; Randel and Lait 1991; 
Manney 1991; Fraser et al. 1993; Lawrence et al. 1995; 
Lawrence and Randel 1996) and theoretical (Hartmann 
1983; Manney et al. 1988; Manney 1991; Manney and 
Randel 1993: Bowman and Chen 1994; Orsolini and 
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Simon 1995) studies have examined the 4-day phenom- 
enon. 

The wave has been observed in both hemispheres but 
is stronger in the Southern Hemisphere (Venne and Stan- 
ford 1982). The horizontal phase structure is variable, 
with both poleward and equatorward momentum flux 
events (Venne and Stanford 1982; Prata 1984; Manney 
1991). The wave is generally barotropic, although some 
episodes show strong equatorward heat flux (Randel and 
Lait 1991). The episodes studied by Randel and Lait 
(1991) and Lawrence and Randel ( 1 996) reveal a dipole 
structure of positive and negative Eliassen-Palm flux 
divergence associated with the 4-day wave, with the 
positive region overlapping negative quasigeostrophic 
potential vorticity gradient, a signature of instability 
processes. 

Theoretical studies have focused on instability of the 
jet structure as the mechanism for 4-day wave growth. 
Hartmann (1983) and Manney et al. (1988) showed that 
periods of 4-day wave growth are consistent with bar- 
otropic instability of the stratospheric polar night jet. 
Manney and Randel (1993) used a 3D instability model 
to show that realistic 4-day wave growth rates exist in 
monthly mean background wind states only if both bar- 
otropic and baroclinic components are included. Their 
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model also predicted a double-lobed structure in the 
perturbation temperature field with one lobe in the upper 
stratosphere and another in the lower mesosphere, with 
rapid phase variation in between. 

The present paper takes advantage of the coverage 
and vertical resolution of the Microwave Limb Sounder 
(MLS) onboard the Upper Atmosphere Research Sat- 
ellite ( UARS ) to examine the 4-day wave during the late 
austral winters of 1992 and 1993. As will be shown 
below, strong 4-day signals exist in MLS temperature, 
geopotential height, and ozone. The vertical coverage 
provided by MLS allows resolution of a double-peaked 
structure in perturbation temperature with an out-of- 
phase relationship between the upper and lower peaks, 
as predicted by Manney and Randel (1993). Spectral 
analyses of MLS ozone reveal a 4-day signal in zonal 
wavenumber 1 near the polar stratopause. A linear ad- 
vective-photochemical model is used to calculate the 
ozone response to temperature and wind perturbations 
associated with the 4-day wave. The vertical variation 
of model ozone amplitude and phase is in good agree- 
ment with observations of MLS ozone. Negative regions 
of quasigeostrophic potential vorticity (PV) gradient and 
positive Eliassen-Palm flux divergence are shown to 
occur, consistent with wave forcing by instability dy- 
namics. 

The 4-day feature is also shown to exist in quasi- 
geostrophic PV derived from MLS geopotential heights. 
The three-dimensional structure of this signal resembles 
the PV “charge” discussed by Hoskins et al. (1985) and 
more recently by Bishop and Thorpe (1994). The latter 
authors have also extended the concept to Ertel PV 
(Thorpe and Bishop 1995), although the present paper 
focuses on quasigeostrophic PV. In this electrostatics 
analogy, a PV anomaly in the atmosphere acts like an 
electrical charge in a dielectric material, with an asso- 
ciated vector field (independent of static stability, den- 
sity, and boundary conditions) that produces “action at 
a distance.” Isolated PV charges induce geopotential 
anomalies, circulation about the vertical axis, and a ver- 
tically oriented temperature dipole similar to observa- 
tions of the 4-day wave presented in this study. 

2. Microwave Limb Sounder data 

a. The MLS instrument 

The Microwave Limb Sounder is 1 of 10 instruments 

on board the Upper Atmosphere Research Satellite . 

MLS uses a 1.6-m scanning antenna to observe the at- 
mospheric limb emission simultaneously in spectral 
bands at 63, 183, and 205 GHz (Barath et al. 1993). 
These observations allow the determination of chlorine 
monoxide, ozone, water vapor, nitric acid, sulfur di- 
oxide, temperature, and geopotential height in the strat- 
osphere and mesosphere during both day and night and 
in the presence of stratospheric clouds and aerosols. This 
study analyzes temperature, geopotential height, and 


ozone data from two austral late-winter time periods, 
14 August-20 September 1992 and 9 August- 16 Sep- 
tember 1993, when MLS was viewing southward (cov- 
erage from 80°S to 34°N). 

b . Temperature data 

MLS temperature is obtained by viewing limb emis- 
sion of 0 2 at 63 GHz. The version 3 data (used here) 
are retrieved on every other UARS pressure level n (the 
UARS level-3 pressure grid is 1000 X 10 "' 2 * * * 6 hPa, where 
n = 0,1,2,- • •) and have precisions ranging from 1.4 K 
at 22 hPa to 3.1 K at 0.46 hPa. MLS temperatures are 
used at eight pressure levels from 22 hPa to 0.1 hPa 
(22, 10, 4.6, 2.2, 1.0, 0.46, 0.22, 0.1 hPa). The last two 
levels are contaminated by large systematic errors and 
are biased by a zonal mean climatology. Wave structures 
estimated at these levels are attenuated and should be 
viewed with some caution. Validation information on 
version 3 MLS temperature is available in Fishbein et 
al. (1996). 

c. Ozone data 

MLS ozone retrievals are obtained from limb emis- 
sion in two bands centered near 205 and 183 GHz. The 
version 3 ozone data are retrieved on every other UARS 
pressure level. The 205-GHz data have smaller uncer- 
tainties throughout the stratosphere, while the 1 83-GHz 
data are less noise limited in the mesosphere. Both 205- 
and 1 83-GHz ozone data are used in this paper from 22 
to 0.22 hPa. Validation information on version 3 MLS 
ozone is available in Froidevaux et al. (1996) and Ri- 
caud et al. (1996). Single-profile precision for ozone 
retrievals is in the 0.1 -0.4 ppmv range, depending on 
altitude. 

d. Geopotential height data 

The MLS geopotential height dataset used in this 
study is a preliminary dataset, which is not available in 
version 3 but is planned in future versions. The geo- 
potential height is obtained by measuring pressure as a 
function of tangent-point height. Pressure is obtained 
from the same radiances as temperature, but the radi- 
ances are sensitive to pressure from 10 hPa to 0.046 
hPa. Tangent-point altitude is determined from the ori- 
entation of the MLS antenna, measured by an angle 
encoder, the UARS attitude control system, and the 
ground-tracked UARS orbit determination. The pressure 
(in log pressure coordinates with atmospheric scale 
height 7 km) precision ranges from 5 m at 10 hPa to 
10 m at 0.1 hPa, while the altitude precision is altitude 
independent and is estimated to be around 100 m. Errors 
in geopotential height (excluding precision) arising from 
the tangent-point altitude are synchronized to the UARS 
orbit and contaminate the zonal mean and high fre- 
quencies near the diurnal tide. However, the analyses 
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Fig. 1 . (a) East-west power spectrum for wave 1 MLS temperature (T nlneved ) at 72°S for 18 August-6 September 
1992. Dashed contours are at 10, 20, and 30 K 2 day, and solid contours are from 40 to 100 K 2 day by 10 K : day. 
(b) Same as (a) but for temperature derived from MLS geopotential height (c),(d) Same as (a),(b) but for 

2-16 September 1993. 


performed in this paper (zonal wave 1, 3-5-day periods) 
have sufficiently long periods to avoid significant effects 
from these errors. Good agreement is found on these 
spatial and temporal scales between MLS temperature 
(coined T ntricved in this paper) and temperature derived 
from MLS geopotential height data (T gcohl ), as discussed 
in section 4a. 

3. Analysis procedure 

MLS data used in this analysis were produced by the 

asynoptic mapping method (Salby 1982a, b; Lait and 
Stanford 1988a). The specific technique used for this 
study is presented in detail in Elson and Froidevaux 
(1993). Briefly, data are binned by latitude and separated 
into ascending and descending time series. The time and 
space axes are rotated to a new coordinate system where 
a fast Fourier transform (FFT) is applied to calculate 
the spectral coefficients. An inverse transformation and 


interpolation allows data to be mapped onto synoptic 
maps, which are processed every 12 hours. Satellite 
orbital period drift generally limits the procedure to ap- 
proximately 1-week analysis periods. The procedure 
here is essentially a remapping in time so that periods 
longer than a week can be analyzed. Data are interpo- 
lated onto horizontal grids of 4° latitude by 5° longitude 
and onto the standard UARS pressure levels. 

The twice-daily synoptic maps described above were 
decomposed in ace and time by applying a two-di- 
mensional sinusoidal transform [see appendix of Ziemke 
and Stanford (1990)]. Power spectral density is defined 
in this study as (2A/)~'[A 2 (fc, co) + B 2 (k, c o )], and am- 
plitude is defined as [A 2 (k, ay) + B 2 (k, a>)] 05 , where A / 
is the unit bandwidth (NDAT A f)" ! . NDAT is the number 
of data points in the time series (NDAT=38 for 1992, 
NDAT =30 for 1993) and A t is the temporal sampling 
interval (0.5 day); A(k, co) and B(k, (o) are the cosine and 
sine coefficients of the discrete Fourier expansions 
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Fig. 2. Time vs longitude plot of T re(rjeved at 64°S and 1 hPa {48 
km) filtered for wave 1 and 3-5 day periods (eastward and westward) 
for 14 August-20 September 1992. Dashed contours are from -4 to 
0 K by I K. Solid contours are from 1 to 4 K by 1 K. 
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Here T, <J>, and fi are temperature, geopotential, and 
ozone mixing ratio; +(-) refers to westward (eastward) 
components; t is time; x is zonal distance; Ax is the 
zonal grid interval of 27ra cos0/72; k is the zonal wave- 
number; oj is the angular frequency; a is the earth’s 
radius; and </> is latitude. A 1-2-1 running mean was 
applied in frequency space to power spectral density, 
but not to amplitude. The time mean was removed be- 
fore applying the FFT to focus the analysis on traveling 
waves. 

Filtered time series in this study were produced by 
first applying the bandpass filter response function (with 
half-amplitude points at 3 and 5 day periods) given by 


Murakami (1979) to the temporal FFT coefficients and 
then recombining the coefficients through an inverse 
FFT. 

As a comparison to the MLS temperature measure- 
ments, temperature was also calculated from MLS geo- 
potential height using the hydrostatic approximation T 
- R~'d< E>/d(ln p), with R = 287 J kg' 1 K 1 (dry air gas 
constant) and constant vertical spacing d ( In p) = 0.387 
(determined from UARS standard pressure surfaces). 
Quasigeostrophic potential vorticity q ' and potential 
vorticity gradient q~ v were calculated using the spherical 
formulation of Matsuno (1970): 
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Here u is the zonal wind, p 0 = p^e~ :JH , p^ is the density 
at the surface (1000 hPa), z is the log-pressure coor- 
dinate, A is longitude, / is the Coriolis parameter, and 
fl is the earth’s angular frequency. Overbars denote zon- 
al means, primes denote deviation from the zonal mean, 
and subscripts denote partial derivatives. The Brunt- 
Vaisalla parameter and atmospheric scale height were 
given the constant values N 2 = 4 X 10 -4 s" 2 and H = 
7 km; q ( was calculated from MLS geopotential height 
data, while q v was calculated from U. K. Meteorological 
Office (UKMO) zonal winds. MLS data were not used 
in the calculation of q y due to the known inaccuracy in 
the zonal mean geopotential height (see section 2d). 

UKMO horizontal winds and temperatures on UARS 
pressure surfaces, obtained from the Goddard Space 
Flight Center Distributed Active Archive Center 
(DAAC), were interpolated with a cubic spline onto 
horizontal grids of 4° latitude by 5° longitude. These 
winds and temperatures were used in the calculation of 
q x and also in the tracer model discussed in section 4c. 

The quasigeostrophic Eliassen-Palm (EP) flux vector 
and the EP flux divergence were calculated with the 
following spherical formulation: 

F = (/> F.) = Po a cos JjfijjT) (4) 


V F = — — -^-(cos^F*) + (5) 

a cos <pd<p dz 

Here u f and v' (meridional wind) are calculated from 
MLS geopotential heights (using the geostrophic ap- 
proximation) and V is from MLS-retrieved tempera- 
tures. In the EP flux diagram (Fig. 1 1), the vector com- 
ponents are multiplied by 2m cos <f>. This factor ac- 
counts for the spherical geometry of the earth to make 
the arrow pattern look nondivergent if and only if V F 
— 0 (see Edmon et al. 1980). 
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4. Results 

a. Temperature 

Space-time spectral decomposition is a useful tool 
for isolating the 4-day wave from other atmospheric 
signals. Figure la presents a periodogram for zonal 
wavenumber 1 power spectral density (defined in sec- 
tion 3 and hereafter called power) for MLS-retrieved 
temperature (hereafter 7 retrievcd ) at 72°S latitude for 18 
August-6 September 1992. A clear 4-day eastward pe- 
riod signal exists near 1 hPa (48 km), and a somewhat 
weaker signal occurs near 0.22 hPa (60 km). These sig- 
nals are also evident in temperature derived from MLS 
geopotential height (hereafter 7 gcohl ) in Fig. lb, although 
the lower signal maximizes near 2.2 hPa (42 km). An 
even stronger 4-day signal is found during 2-16 Sep- 
tember 1993 in both r retrieved (Fig. lc) and 7 gcohl (Fig. 
Id). The 4-day signal dominates the traveling wave 
spectra for this time period with two lobes centered near 
2.2 hPa (42 km) and 0.22 hPa (60 km). Sections 4d and 
4f further discuss the validity of the 4-day signal in 
terms of uncertainty and statistical significance. 

The 1992 4-day wave event was further isolated by 
filtering a 38-day series of 7 relrieved for zonal wavenumber 
1 and periods (both eastward and westward) of 3-5 days 
and plotting on a Hovmdller (time vs longitude) diagram 
for 64°S and 1 hPa (Fig. 2). The 4-day wave amplitude 
increases over the first 10 days and remains constant 
for another 10 days before dissipating near 1 September 


1992. The mode makes at least five revolutions around 
the pole with a rotation period near 3.8 days. 

Spectral plots of T remevcd for zonal wavenumbers 2 and 

3 (not shown here) for the 1992 event show signals at 
1 hPa with eastward periods of 2.0 and 1.3 days, re- 
spectively. These sinusoidal components move at the 
same phase speed of the wave 1 , 4-day signal and com- 
bine with wave 1 to form a “warm pool” of air at 1 
hPa (coupled with a slightly less prominent “cold 
pool”) that rotates around the pole with a period near 

4 days. This feature is shown in Fig. 3, where 7 retricved 
data have been filtered over 14 August-20 September 
1992 for waves 1-3 and eastward rotation periods of 
3-5 days. (Wave 4 and higher were not analyzed because 
the frequency of interest, 4 days divided by integer 
wavenumber, is near or outside the Nyquist limit of 
approximately 1 day 1 .) The warm pool rotates around 
the pole as an identifiable entity throughout the record 
shown displaying the quasi-nondispersive (soliton-like) 
nature of the feature. 

The perturbation temperature (7 retneved ) amplitude for 
wave 1 and 3.8-day eastward period is plotted in Fig. 
4a for the 1992 4-day wave event. A large lobe exists 
in the upper stratosphere, peaking near 48 km and 72°S 
with maximum amplitude of 4.2 K. A smaller upper 
lobe is centered at 60 km and 68°S. This double-lobed 
structure is similar to that predicted by the instability 
study of Manney and Randel (1993). In their results, 
the upper lobe is larger and extends well into the mes- 
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(a) Tretrieved (1992) 
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Fig. 4. (a) F remcvcd amplitude for wave 1 and 3.8-day eastward period 
for 18 August-6 September 1992. Contours are every 10% interval 
of the maximum (4.2 K), with values higher than 50% shaded, (b) 
Tranewd amplitude for wave 1 and 3.75-day eastward period for 2-16 
September 1993. Contours are every 10% of the maximum (5.8 K), 
with values higher than 50% shaded. 


osphere. However, r retrieved is largely climatology at 66 
km, forcing zonally asymmetric wave amplitudes to 
zero; therefore, the abrupt fall-off with height above 60 
km in Fig. 4a may not be real. The 1993 series also 
displays a double-lobed structure (Fig. 4b), with maxima 
near 42 km, 76°S (amplitude of 5.8 K) and 60 km, 80°S. 

Manney and Randel (1993) also predict rapid vertical 
phase variation between the two lobes. Such behavior 
can be explicitly seen in height versus longitude plots 
of r relneved at 72°S, filtered for waves 1-3 and eastward 
rotation periods of 3-5 days. Figure 5a presents 8 days 
during the peak amplitude of the 1992 time series. The 
plots reveal a warm anomaly near 48 km bounded above 
by a cold anomaly near 60 km. This dipole pattern con- 
tinues throughout the 4-day wave life cycle, suggesting 


that the two peaks constitute one entity. The 1993 series 
also reveals this vertical temperature dipole (Fig. 5b). 
Here the lower peak has a strong wave 1 structure that 
begins near 48 km and descends to 42 km over this 
8-day period. The upper feature, which maximizes near 
60 km, is approximately 180° out of phase with the 
lower feature. The quadrupole pattern produced by the 
two out-of-phase lobes remains stable for at least two 
weeks, through the end of the UARS month (16 Sep- 
tember 1 993), when MLS changed to northward viewing 
(coverage from 80°N to 34°S). 

b. Geopotential 

The vertical temperature dipole structure seen in Fig. 
5 suggests that the perturbation geopotential amplitude 
should increase with height throughout the region of the 
lower temperature lobe and decrease with height in the 
region of the upper lobe. This is confirmed in Fig. 6a, 
which shows the geopotential amplitude for the 1992 
series for eastward moving wave 1 with a period of 3.8 
days. The signal maximizes near 57 km, at the boundary 
between the warm and cold anomalies (see Fig. 5a), and 
decreases above it. The 1993 geopotential amplitude 
(Fig. 6b) peaks near 48 km, consistent with the pertur- 
bation temperature structure for that series (Fig. 5b). 
The geopotential signal decreases with altitude more 
rapidly in the 1993 time series, due to a larger- amplitude 
upper temperature lobe in 1993 than in 1992. The me- 
ridional structure shown here is similar to previous ob- 
servations of the 4-day wave in geopotential height 
(Manney 1991; Lawrence and Randel 1996). 

c . Ozone 

Other than the unpublished conference report by Ran- 
del et al. (1992), the present paper is the first to present 
observations of the 4-day wave in ozone. Periodograms 
for wave 1 MLS ozone (both 205- and 183-GHz data) 
at 72°S during 18 August-6 September 1992 are pre- 
sented in Fig. 7. A 4-day peak exists near the stratopause 
(1.0 hPa), coincident with the r rtflricved maximum at the 
same altitude, see Fig. la. A 2-day eastward period sig- 
nal exists at 2.2 hPa; the source of this signal is unknown 
and is not thought to be related to the 4-day wave. The 
2-day feature does not exist in spectra from 5 to 20 
September 1992. 

To ascertain the relationship between the 4-day fea- 
ture in ozone and the other atmospheric variables, a 
linear advective-photochemical model was used to cal- 
culate the ozone response to temperature and meridional 
wind perturbations associated with the 4-day wave. The 
procedure combines the linearized thermodynamic and 
ozone continuity equations to obtain an expression for 
the ozone mixing ratio (Hartmann and Garcia 1979; 
Rood and Douglass 1985; Randel 1990, 1993). Because 
the 4-day wave amplitudes in the various components 
(T, u, v, and p) are observed to be much less than their 
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(a) Tretrieved W 1-3, 3-5 day (22-29 August 1992) 
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Fig. 5. (a) 7' relnevcd on a height vs longitude grid at 72°S filtered for waves 1-3 and 3-5-day eastward rotation period for 22-29 August 
1992. (b) Same as (a) but for 4-11 September 1993. Contours are every 1 K with solid (dashed) line having positive (negative) values. 
Here, W (C) refers to warm (cold) anomalies. 


respective zonal mean values {T C F, etc.), we use a 
linearized thermodynamic equation, 

T ; + uT x + v f T y + w'S = 0 

S = j^IkT + Ht z l (6) 

where w r is the vertical wind component, k = Rlc p ~ 
2/7, and H = 7 km; and a linearized ozone continuity 
equation, 

fi\ + uyJ x + v + w'/I, = -T^ r - 0f. (7) 


The photochemical relaxation rate, F, and the linearized 
ozone response to small temperature perturbations, ©, 
are shown in Fig. 8a, adapted from Hartmann and Garcia 
(1979), who studied ozone transport by planetary waves 
in the winter stratosphere using a linearized /3-plane 
model centered at 60°. Radiative damping of the tem- 
perature perturbations is neglected, since the radiative 
timescale is of order 10 days in the region of interest 
(26-57 km) [see Hartmann and Garcia (1979)], while 
the advective timescale [k(u — c)]" 1 is of order 1-2 days 
from 26 to 42 km, and the photochemical timescale T -1 
is less than 1 day from 42 to 57 km. Radiative heating 
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(a) Geopotential (1992) 
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(b) Geopotential (1993) 
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Fig. 6. (a) MLS geopotential amplitude for wave 1 and 3.8-day 
eastward period for 18 August-6 September 1992. (b) Same as (a) 
but for 2-16 September 1993, 3.75-day eastward period. Contours 
are at every 10% interval of the maximum with values higher than 
50% shaded. 


from ozone perturbations has also been neglected, as in 
Randel (1990, 1993), since it is likely small compared 
to changes due to advection. Solving Eq.(6) for w' and 
substituting into Eq.(7) gives the following differential 
equation for fx!\ 


^ + u a~ + r )a + 

dt dx 


Mv 



[S\dt 



_ 0 


r 


= o. 


( 8 ) 


To relate p! to T and v’ , we write (/x\ T\ v ) - (M, 
T, V) e ikix " v, \ where M, T, and V denote complex-valued 
perturbation amplitudes. From this, 


M = —-r—— + — (9, 

T + ia 

£ = fi y - fiJJS a = k{u - c) b = k{u - c)fiJS. 

Figure 8b presents the resulting model ozone ampli- 
tude for the wave 1, 3.8-day period signal at 60°S (thick 
solid line). Also plotted are the observed ozone ampli- 
tudes at 60°S for 205-GHz (thin solid) and 183-GHz 
(dashed) data. Error bars on the observed data were 
calculated using the procedure described in section 4d. 
Reasonable agreement is found between the observed 
and model amplitudes, with low- values in the middle 
stratosphere and lower mesosphere, and a maximum in 
the upper stratosphere (near 42 km). Results from a 
purely advective model (T = 0 = 0) and a purely 
photochemical model [fi amp = (0/T) T amp , see Randel 
( 1 990)] show a transition region between dynamical and 
photochemical control from approximately 37 to 49 km, 
where both dynamics and photochemistry play an im- 
portant role. 

The phase of the model and observed 4-day feature 
is provided in Fig. 8c. Good agreement exists between 
the 205-GHz observed phase and the model throughout 
the stratosphere, while the 183-GHz observed phase 
agrees with the model in the lower mesosphere but di- 
verges somewhat in the stratosphere. As mentioned in 
section 2c, the 205-GHz data have smaller uncertainties 
in the stratosphere, while the 183-GHz data are less 
noise limited in the mesosphere, so the phase results 
showm here are consistent with the data quality. 

The good agreement (at least to first order) between 
model and observations suggests that the observed 
4-day signal in ozone is produced largely by the coupled 
effects of constituent advection by the 4-day wave per- 
turbation winds in the presence of zonal-mean constit- 
uent gradients and photochemical processes associated 
with the 4-day wave perturbation temperature and mix- 
ing ratio anomalies. 

d. Uncertainty in 4-day wave amplitude 

In order to assess the uncertainty in the 4-day wave 
amplitude calculation, a random number generator was 
used to make a noise dataset with level 3AT format 
(along-track, evenly spaced in time) with values be- 
tween -0.5 and +0.5 (arbitrary units), giving a root 
mean square (rms) variation of 0.29. This represents the 
random noise in a given MLS measured quantity with 
an rms error (defined here as uncertainty) of 0.29. These 
data were passed through the asynoptic mapping routine 
described in section 2 and synoptic maps were produced. 
The rms variation in the synoptic maps is 0.19, which 
is 65% of the input value. The decrease is mainly due 
to interpolation to a regular horizontal grid. Fourier de- 
composition was then applied to isolate the amplitude 
of the wave 1, 3.8-day component from 500 different 
time series, each randomly generated. This Monte Carlo 
procedure yielded rms variations in amplitude of 0.018 
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(a) 205 GHz Ozone (b) 183 GHz Ozone 

Perfod [days] Period [days] 




Fig. 7. (a) East-west power spectrum for wave 1 MLS ozone (205 GHz) at 72°S for 18 August-6 September 
1992. Contour increment is 0.02 ppmv-day. Dashed lines are from 0.02 to 0.08 ppmv- day. Solid lines are from 
0.10 to 0.20 ppmv 2 day. (b) Same as (a) but for 183-GHz MLS ozone. 


for the 1992 simulated series (38 time components used, 
72 longitudes) and 0.020 for the 1993 simulated series 
(30 time components used, 72 longitudes). The uncer- 
tainty in amplitude is 0.018/0.29 - 1/16 times the un- 
certainty in the level 3 AT data for 1992 and 1/14 for 
1993. The difference between 1/16 and 1/14 is due to 
the different bandwidth used in each year. This method 
can now be used to obtain estimated uncertainties in 
4-day wave amplitudes observed in MLS temperature 
and ozone. 

The estimated uncertainty in version 3 MLS temper- 
ature at 1 and 2.2 hPa is 3.4 and 2.2 K, respectively 
(Fishbein et al. 1996). The estimated wave 1, 3.8-day 
eastward period amplitude uncertainty in the lower peak 
of Fig. 4a (1992) is 3.4/16 = 0.21 K, which is much 
smaller than the observed amplitude of 4.2 K. In 1993 
the uncertainty in the lower peak (Fig. 4b) is 2.2/14 = 
0.16 K, which is much less than the observed amplitude 
of 5.8 K. This Monte Carlo analysis shows that the 
observed amplitudes are not significantly affected by 
random noise in the data. 

This method was also used to obtain estimated errors 
for 205- and 183-GHz MLS ozone data. The error bars 
shown in Fig. 8b were computed using the above for- 
mulation along with data uncertainties at different pres- 
sure levels as given by Froidevaux et al. (1996). 

e. Statistical assessments 

Because the 4-day wave is a well-established atmo- 
spheric phenomenon (see references in the introduc- 
tion), statistical assessment of the 4-day features found 
here can be made on the basis of a priori statistics (Mad- 
den and Julian 1971). Briefly, before analysis it is de- 
cided to examine spectral features near 4-day periods. 
[If spectral peaks are examined at previously unreported 


frequencies, a posteriori statistics must be used (Madden 
and Julian 1971).] The null hypothesis is made that any 
spectral peak at 4 days is due to random sampling fluc- 
tuations. Under this hypothesis, the ratio of spectral peak 
to nearby spectral background can be estimated statis- 
tically. For the Fourier transform spectra used here, chi- 
squared statistics are employed to estimate the maxi- 
mum peak to background ratio that will occur in 95% 
of the cases for random sampling fluctuations. If the 
observed peak exceeds this ratio, with 95% confidence 
the spectral feature can be said to be statistically sig- 
nificant. 

Figures 9a, b show the T retrieved spectra at 2 and 0.2 
hPa during 2-16 September 1993. These are plotted 
without smoothing in Figs. 9a, b and correspond to the 
smoothed plot shown in Fig. lc. Also indicated by 
dashed lines in Fig. 9 are the estimated background and 
95% confidence limit (for two degrees of freedom, ap- 
propriate to this case of unsmoothed Fourier transform 
spectral estimates). Both features easily satisfy a priori 
statistical significance. 

The 1992 temperature spectral 4-day peak, shown 
smoothed in Fig. la, reveals a strong peak centered at 
1 hPa that also easily exceeds the a priori 95% confi- 
dence limit (not shown here). The weaker temperature 
peak at 0.22 hPa is close to a stronger lower-frequency 
feature. By examination of unsmoothed spectral results 
(not shown here), it can be shown that the 0.22-hPa 
4-day peak also meets the criterion for a priori statistical 
significance. 

The 183-GHz ozone spectrum at 1 hPa for 1992 
(shown smoothed in Fig. 7b) is plotted unsmoothed in 
Fig. 9c. The 4-day peak exceeds the 95% confidence 
limit and is judged statistically significant. The 205-GHz 
ozone 4-day peak at 1 hPa (Fig. 7b) also exceeds the 
95% confidence limit (not shown here). 
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(a) 0 [io‘ ppmv s' 1 K'J (a) Tretrieved (1993) at 2.0 hPa 




(b) Amplitude 


(b) Tretrieved (1993) at 0.2 hPa 




(c) Phase 



Fig. 8. (a) Ozone photochemical parameters adapted from Hart- 
mann and Garcia (1979). The curve Tis the photochemical relaxation 
rate, and 0 is the linearized ozone response to small temperature 
perturbations, (b) Wave l ozone amplitude (in ppmv) at 60°S for 3.8- 
day eastward period from 18 August to 6 September 1992. Lines are 
model (thick solid), 205-GHz observations (thin solid) and 183-GHz 
observations (dashed). Error bars were calculated using the procedure 
described in section 4d. Large (small) caps show the extent of the 
205 (183) GHz ozone error bars, (c) Same as (b) but for phase and 
excluding error bars. 


(c) 183 GHz Ozone (1992) at 1.0 hPa 



Fig. 9. Unsmoothed power spectra for wave 1 at 72°S plotted along 
w ith estimated background and 95% confidence level, (a) T Mnewi for 
the 1993 series at 2 hPa. (b) 7 , r(lrieveU for the 1993 series at 0.2 hPa. 
(c) 183-GHz ozone for the 1992 series at 1.0 hPa. 

/. Zonal wind, q r and EP flux structure 

Theoretical studies show that barotropic and/or bar- 
oclinic instability of the polar night jet can produce 
quasi-nondispersive modes similar to the observed 
4-day wave. The instability episodes are generally ac- 
companied by regions of negative quasigeostrophic po- 
tential vorticity gradient q y and nonzero EP flux diver- 
gence. Figure 10a shows the zonal mean UKMO east- 
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(a) UKMO Zonal Wind 18 August - 6 September 1992 



(b) QGPV Gradient 18 August - 6 September 1992 



Fig. 10. (a) UKMO zonal mean zonal wind (westerly, in m s -1 ) 
averaged from 18 August to 6 September 1992. (b) Quasigeostrophic 
potential vorticity gradient, q y , averaged from 18 August to 6 Sep- 
tember 1992. Contour increment is 1.0 X 10" "m 's 1 with negative 
values dashed. 


west wind averaged over the life cycle of the 1992 4-day 
wave event (18 August-6 September 1992). A strong 
stratospheric westerly jet is centered near 34 km, 60°S. 
The associated wind curvature produces a region of neg- 
ative q y near 50°S at altitudes above 50 km (Fig. 10b). 
Daily wind maps reveal a “double jet” structure near 
the beginning of the time series, which is often present 
during 4-day wave events (Manney 1991; Randel and 
Lait 1991). Manney et al. ( 1 996) showed that daily maps 
of Ertel PV (on isentropic surfaces) calculated from 
UKMO winds sometimes exhibit unrealistic structure in 
the Southern Hemisphere winter middle and upper strat- 
osphere poleward of about 70°S. Our analyses here, 
however, focus on quasigeostrophic PV averaged over 
20 days. Comparison of Fig. 10a with National Centers 
for Environmental Prediction (NCEP, formerly NMC) 
zonal winds (not shown here) suggests good overall 


EP Flux Wave 1 3.8 day 18 August - 6 September 1992 



Fig. 11. Eliassen-Palm flux diagram calculated for the wave 1, 
3.8-day eastward period components of MLS temperature and geo- 
potential height from 18 August to 6 September 1992. Contours are 
DF = V-F/(p q a cos</>) starting at 0.1 and incremented at 0.4 ms 1 
day" 1 with solid (dashed) line indicating positive (negative) values. 
Arrows are the Eliassen-Palm flux vectors. 

agreement. However, q } calculated from the NCEP data 
(which extends only to 54 km) does not show the neg- 
ative region near 50°S for this time-averaged period, 
although negative q y does exist near 50°S, 54 km for 
several of the individual NCEP-derived daily maps dur- 
ing the first half of this time series. (Note that time- 
averaged fields are more smoothed than the wind field 
on any given day, so they often underestimate instabil- 
ity.) 

Linear theory requires that wave growth by barotropic 
instability be accompanied by negative PV gradients. 
[The integral over latitude of q y times a positive definite 
quantity must be zero for wave growth to occur by 
barotropic instability. However, baroclinic instability is 
possible with positive q y , given suitable variation with 
latitude of the vertical gradient of u at the lower bound- 
ary; see Holton (1992).] The presence here of negative 
q Y regions is consistent with the fact that negative q y is 
a necessary, although not sufficient, condition for wave 
growth by barotropic instability. 

Figure 1 1 displays the EP flux diagram for this wave 
episode, giving details about the net effect on the zonal 
mean state of eddy and momentum fluxes due to wave 
events. The EP flux was calculated using MLS tem- 
perature and geopotential height from 18 August to 6 
September 1992, filtered for the wave 1 and 3.8-day 
eastward period. The geostrophic approximation was 
used to calculate the winds from the 4-day wave per- 
turbation geopotential field. The EP vectors are plotted 
over contours of DF = V F / (p 0 a cos</>), which is a 
force per unit mass on the zonal-mean wind (in m s -i 
day -1 ) due to the eddy effects associated with the wave 
1 component of the 4-day wave. Positive (negative) 
regions of DF act to accelerate (decelerate) the zonal 
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(a) Wave 1 PV (1992) (b) Wave 1 PV (1993) 


Per i od [ days 1 

1.0 1,3 2.0 4.0 4,02.01.31.0 


Period [ days ] 

1.01.32.04.0 4.02.01.31.0 




Westward Freq. [1/day] Eastward 


Westward Freq. !l/day) Eastward 


Fig. 12. (a) East-west power spectrum for wave 1 quasigeostrophic potential vorticity derived from MLS 
geopotential height at 72°S for 18 August-6 September 1992. (b) Same as (a) but for 2-16 September 1993. 
Contours are arbitrary with solid lines representing larger values. 


mean westerly wind. Positive DF exists between 50° 
and 75°S at the top of the plot, near 57 km, acting to 
accelerate the lower-mesospheric winds, while a nega- 
tive DF region near 54 km, 72°S acts to decelerate the 
stronger stratopause winds (see Fig. 10a). Figure 11 
predicts eddy-induced reduction of vertical shear in the 
60°-75°S lower-mesosphere region, consistent with bar- 
oclinic instability mechanisms driving the 4-day wave 
in the lower mesosphere. The EP flux structure above 
48 km is similar to the model results from Manney and 
Randel (1993, their Fig. 2d), with a region of negative 
DF bounded above and to the north by positive DF. 

The EP flux vectors indicate strong equatorward heat 
flux (downward pointing arrows) and equatorward mo- 
mentum flux (poleward pointing arrows) above 48 km 
and poleward of 65°S (also similar to Manney and Ran- 
del 1993, Fig. 2d) and largely poleward momentum flux 
elsewhere. The vectors are predominantly vertical (in- 
dicating baroclinic processes) in the lower mesosphere, 
while strong momentum flux (indicated by horizontal 
components of EP flux vectors) near 42 km, 60°S sug- 
gests that barotropic processes play a sizable role there. 
This is confirmed by the positive (70°S) and negative 
(60°S) DF regions at 42 km, which act to reduce the 
latitudinal shear associated with the stratospheric jet (a 
signature of barotropic instability). 

Manney and Randel (1993) and Lawrence and Randel 
( 1 996) show overlapping regions of negative q v and pos- 
itive V-F during 4-day wave episodes. This spatial co- 
incidence is recognized as a source region of wave ac- 
tivity A = (Vt)q n fq y , where q represents the perturba- 
tion quasigeostrophic potential vorticity. Neglecting 
diabatic and frictional effects, a conservation equation 
can be written (Andrews and McIntyre 1976, 1978): 

r)A 

F V-F = 0. (10) 

dt 


Clearly, when q y is negative, the wave activity will in- 
crease in magnitude when V-F > 0. One region of neg- 
ative q y occurs in the vicinity of positive DF near 57 
km, 52°S. This is consistent with instability dynamics 
driving the 4-day wave in the lower mesosphere. Ac- 
cording to Eq. (10), wave activity will also increase if 
q v > 0 and V-F < 0. The regions of negative DF near 
54 km, 80°S-65°S and 42 km, 60°S overlap positive q v% 
suggesting instability forcing may also exist in these 
regions. 

g. Quasigeostrophic potential vorticity 

Spectral analyses of quasigeostrophic PV calculated 
from MLS geopotential heights reveal 4-day signals in 
both the 1992 and 1993 series. Figure 12 shows per- 
iodograms for 72°S during 18 August-6 September 
1992 (Fig. 12a) and 2-16 September 1993 (Fig. 12b). 
Clear 4-day eastward period signals in PV are seen near 
the stratopause (48 km or 1 hPa) during both time pe- 
riods. (In 1992, a significant longer-period feature is 
seen in the upper stratosphere. This eastward moving 
power with periods of several weeks is not directly re- 
lated to the 4-day wave considered here.) A height ver- 
sus latitude plot of PV amplitude for wave 1 and 3.8- 
day eastward period is provided in Fig. 13a for the 1992 
series. A region of large PV amplitude is centered near 
50 km, 68°S. Figure 13b shows the wave 1, 3.75-day 
signal in 1993, which has a similar structure to the 1992 
feature, although centered slightly lower and with con- 
tours somewhat more compact. 

The three-dimensional structure of the perturbation 
PV resembles the “PV charge” concept discussed by 
Hoskins et al. (1985) and more recently by Bishop and 
Thorpe (1994), who considered the effects on the at- 
mosphere of a uniform spherical quasigeostrophic PV 
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(a) Wave 1 PV (1992) 



Latitude [degrees] 


(b) Wave 1 PV (1993) 



Latitude [degrees] 

Fig. 13. (a) Quasigeosirophic potential vorticity amplitude for wave 
1 and 3.8-day eastward period for 18 August-6 September 1992. (b) 
Same as (a) but for 2-16 September 1993, 3.75-day period. Contours 
are at every 10% interval of the maximum with values higher than 
50% shaded. 

anomaly in a basic state with constant static stability 
and density and in a two-layer vertical static stability 
basic state representing the effects of the tropopause. In 
this analogy, the atmosphere represents an isotropic di- 
electric material, static stability represents dielectric 
properties, geopotential (or streamfunction) represents 
electric potential, and an isolated PV anomaly represents 
an electrostatic charge, with an accompanying vector 
field that produces action-at-a-distance (a fundamental 
concept in PV thinking). This “PV field” (which is 
independent of static stability, density, and boundary 
conditions) induces perturbation temperature, geopoten- 
tial, and wind fields as shown in Fig. 14a (for the case 
of constant static stability), adapted from Figs. 1 and 2 
of Bishop and Thorpe (1994). Here a negative spherical 
PV charge induces a vertical temperature dipole with 
cold temperature perturbation on top, spherical geopo- 


(a) PV Charge Schematic 



(b) 1992 4 — Day Wave 



Latitude [degrees] 

Fig. 14. (a) “PV charge” schematic adapted from Bishop and 
Thorpe (1994), Figs. 1 and 2. The negative PV anomaly induces 
temperature, geopotential, and wind (solid line with arrows) anom- 
alies as shown. Note: the PV, temperature, and geopotential anomalies 
are viewed here as a meridional cross section of the three-dimensional 
PV charge, while the circulation is in a horizontal plane that projects 
out of the page, (b) Amplitude of wave 1, 3.8-day eastward period 
signal from 18 August to 6 September 1992. Hatched region includes 
PV amplitude greater than 90% of the maximum. Solid lines are 
perturbation temperature amplitude of 60%-90%, incremented by 
10% of the maximum. Dashed lines are perturbation geopotential 
amplitude at every 10% of the maximum. Adapted from Figs. 4a, 
6a, and 13a of this paper. 


tential anomaly, and anticyclonic circulation about the 
vertical axis. Note that circulation shown here is for the 
Southern Hemisphere (f < 0); the direction should be 
reversed for the Northern Hemisphere (f > 0). 

In comparison, Fig. 14b shows the amplitude of the 
wave 1, 3.8-day period mode during the 1992 series as 
a function of height and latitude. The hatched region 
includes PV amplitude greater than 90% of the maxi- 
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mum. Solid lines are perturbation temperature amplitude 
of 60%-90% of the maximum, incremented by 10%, 
and the dashed lines are perturbation geopotential am- 
plitude at every 10% of the maximum. Here we see the 
PV anomaly bounded above and below by two tem- 
perature lobes, which are out of phase, resembling the 
cold and warm perturbations in Fig. 14a. The geopo- 
tential anomaly surrounds the region of large PV am- 
plitude, although centered slightly higher. There would 
also be an induced circulation about the vertical axis 
associated with the geopotential anomaly, similar to that 
shown in Fig. 14a. 

This analogy suggests that the 4-day wave may be 
more fundamentally viewed as a PV anomaly than a 
temperature anomaly. The wave life cycle begins with 
an unstable atmosphere that allows growth in PV anom- 
aly amplitude, which in this paradigm induces anomalies 
in temperature, geopotential. and wind. This view also 
appeals to the “wave activity’' interpretation of eddy 
features discussed in section 4f, since PV anomalies 
relate directly to the EP flux divergence forcing of the 
zonal-mean state. The PV charge paradigm, coupled 
with its generation by instability mechanisms, appears 
to provide a satisfying physical explanation for the 
4-day wave. 

5. Summary 

This study takes advantage of the vertical coverage 
and resolution of the Microwave Limb Sounder instru- 
ment to analyze 4-day wave features during August- 
September 1992 and 1993. Strong 4-day wave events 
are found during both time series. A double-peaked 
structure is found in temperature with maxima in the 
upper stratosphere and lower mesosphere and strong 
vertical phase variation in between, as predicted by the 
model of Manney and Randel (1993). Four-day eastw ard 
period signals are also evident in MLS ozone, geopo- 
tential height, and quasigeostrophic potential vorticity 
derived from MLS geopotential height. The ozone sig- 
nal compares well with results from a linear advective- 
photochemical model that calculates the ozone response 
to temperature and wind perturbations associated with 
the 4-day wave. A negative region of quasigeostrophic 
potential vorticity gradient coexists with positive EP 
flux divergence, suggesting that instability dynamics 
plays a role in the w'ave forcing. The three-dimensional 
wave structure is shown to resemble the potential vor- 
ticity charge concept, which includes a vertical tem- 
perature dipole and anomalies in PV, geopotential, and 
horizontal wind. The PV paradigm, in which tempera- 
ture, geopotential, and wind anomalies are induced by 
a potential vorticity charge, itself the result of instability 
growth processes, offers an elegant and conceptually 
satisfying explanation of the physics of the 4-day wave. 
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